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Spin incommensurability (IC) has been recently experimentally discovered in the hole-doped 
Ni-oxide chain compound Y2- x Ca x BaNi05 (G. Xu et al, Science 289, 419 (2000)). Here a two 
orbital model for this material is studied using computational techniques. Spin IC is observed in a 
wide range of densities and couplings. The phenomenon originates in antiferromagnetic correlations 
"across holes" dynamically generated to improve hole movement, as it occurs in the one-dimensional 
Hubbard model and in recent studies of the two-dimensional extended t-J model. The close proximity 
of ferromagnetic and phase-separated states in parameter space are also discussed. 

PACS numbers: 71.10.-w, 75.10.-b, 75.30.Kz 



One of the most remarkable results in hole-doped two- 
dimensional (2D) high temperature superconductors is 
the existence in neutron-scattering experiments of spin 
incommensurability (IC)u This result is compatible with 
the formation of stripes, with a 7r-shift in the antifer- 
romagnetic (AF) order across the stripes. Consider- 
able theoretical work has been devoted to the expla- 
nation of these structures, but no consensus has been 
reached. Recently, a phenomenon analogous to the spin 
IC in cuprates has been reported La the one-dimensional 
(ID) compound Y2_ x Ca x BaNi05.B The parent material 
(x=0) is a spin-one (S=l) chain described by a Heisen- 
berg model, with a spin liquid ground-state. Upon dop- 
ing, the resistivity is drastically reduced,cl and the mag- 
netic IC arises.13 These novel results raise the interesting 
possibility of a common origin of the spin IC phenomena 
observed in doped 2D cuprates and ID nickelates. 

The purpose of this paper is to present theoretical cal- 
culations searching for spin IC in models for doped S—l 
chains. Theoretical studies of these systems have been 
presented before, ErO but the magnetic order upon doping 
has not been investigated in detail. The main result of 
our effort is that spin IC can indeed be observed when 
holes are doped into a spin-gapped S=l chain. The ef- 
fect originates in the local existence of dynamically in- 
duced robust AF correlations between the spins located 
at both sides of a mobile hole ( "across-the-hole" ) . This 
spin structure is generated to facilitate the hole move- 
ment, namely to improve the kinetic energy portion of 
the Hamiltonian. This is qualitatively similar to re- 
sults reported in ID Hubbard models with spin-charge 
separation,!!! and also in studies by Martina,et al. on lad- 
ders and 2D models with very mobile holes. □ Based on the 
results reported here, it is natural to conjecture that the 
spin IC in ID Ni-oxides may originate in similar tenden- 
cies of these systems toward spin-charge separation. In 



the ID S=l/2 system this separation is complete, while 
in 2D S=l/2 and ID S—l spin-gapped systems it ap- 
pears to occur only at short distances. However, this is 
sufficient to induce a robust spin IC. 

The proposed Hamiltonian for doped S=l chains is 
H = - 22 tab(cl aa c jba + h.c.) + U ^ n i^n iai 

(ij)aba ia 
+ U' n ila n l2cr ' - J ^ C- lCT Cji CT 'c|2a' C »2cr, (1) 



where the notation is standard, namely rii aa is the num- 
ber operator at site i, orbital a, and spin cr, a = 1(2) 
corresponds to orbital 3d(3z 2 — r 2 ) (3d(x 2 — y 2 )) of the 
Ni +2 ion, and the chain direction is taken to be the pr- 
axis . The hopping amplitudes between nearest-neighbor 
sites are tu = 4/3, t 22 = h 2 = t 2l = 0.E3 U and U' 
are onsite intra-orbital and inter-orbital Coulomb inter- 
actions, respectively, while J(> 0) is the Hund coupling 
which favors parallel spins on the same site. N ote . t hat 
orbital rotational invariance requires U =U' + J.tlrll3 For 
simplicity, in this first attempt to understand the spin IC 
effect the oxygen ions are not explicitly included in the 
Hamiltonian. The use of Zhang- Rice "doublets!!, (instead 
of singlets) in previous literature of Ni-oxidesQO justifies 
in part this approximation. 

The many-body technique used here to ana- 
lyze ground-state properties of Hamiltonian (Q) is 
the Deusity-Matrix Renormalization Group (DMRG) 
method.Ej The finite-system variation of DMRG was ap- 
plied, working with open boundary conditions (OBC). 
Truncation errors were kept around 10 -5 or smaller, us- 
ing about 100 states. In order to characterize ground- 
state properties, a variety of observables have been mea- 
sured, such as the spin structure-factor defined as 
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where Sj a = J2 a (3 c ]aa (T af3Cjai3, ( ) denotes ground-state 
expectation value, and L is the length of the chain. 
The charge structure-factor was also measured, but no 
evidence of charge ordering has been observed. Also 
the local density, (n*) = ^2 aa {niacr), and total-spin z- 
component, (S*) = 2„(S| ), were studied, but the main 
results were obtained focusing on spin correlations. 
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FIG. 1. Results at n—2 for a 16-site chain, (a) Phase di- 
agram, (b) Spin correlations at (U' , J) = (6, 4). For compari- 
son, the correlations of the AF S=l Heisenberg chain are also 
shown. The inset shows the same but for ({/', J)=(16, 13). 
Details can be found in the text. 

Let us consider density n—2 first, expected to represent 
the undoped material. Fig. la shows the ([/', J) phase 
diagram for a system of 16 sites and 32 electrons. Each 
of the three marked regions have distinct features. In the 
regime above line J—U' the Hund coupling J provides an 
electronic attractive interaction which competes with the 
U' repulsion. In order to have the expected overall repiJ-, 
sive total inter-orbital interaction, J < U' is required.til 
In addition, to create a spin-1 chain, as observed in ex- 
periments, a robust Hund coupling is needed. In fact, 
inside the region labeled S=l it was observed that the 
total spin at site j, defined as Sj =_Sji + Sj2, is to an ex- 
cellent approximation equal to onelij for any site j. This 
is the region of interest in the present problem. In the 
regime labeled S^l, the mean spin is less than 1 since 
J is not sufficiently strong.EZl The frontier S=l-S^l has 
error bars caused by finite size effects, estimated from 
chains with L=8, 16 and 20 sites. 

In the region of local S=l states it is necessary to ver- 
ify whether a Heisenberg spin-1 chain provides a proper 
effective model for the system. For this purpose, ground- 
state spin correlations were calculated for several cou- 
plings inside the S—l region. As shown in Fig. lb, these 
correlations behave quite similarly for two typical sets of 
couplings (U' , J), and in addition they are nearly identi- 



cal to those arising from the much simpler spin-1 Heisen- 
berg chain model. In both cases, the results at a given 
distance were averaged over the entire chain to improve 
the convergence to the thermodynamic limit when differ- 
ent OBC system sizes are used. 




FIG. 2. The spin structure factor for systems of 16 and 20 
sites (indicated in parenthesis) and different electron densi- 
ties: (a) U'=6 and J=4; (b) [7'=16 and J=13. 

Inside the S=l region, it can be shown that not 
only the ground-state properties, but also the low-energy 
physics of Eq.(l) reduces to the spin-1 Heisenberg chain, 
with an effective antiferromagnetic coupling Jneis as an 
overall scale, which depends on U' and J. This interest- 
ing result confirms that our model at n=2 is realistic. In 
fact, for several values of (V, J) and several lattice sizes, 
the effective Heisenberg coupling Jneis was estimated by 
comparing energy gaps of Eq.(l) with those of the spin-1 
Heisenberg chain. For given values of (U 1 , J), data com- 
ing from systems with sizes between 4 and 20 sites yield 
approximately the same value for Jueis ■ Even though a 
closed form for the dependence of Jneis with couplings 
is unknown, it is clear from our results that Jueis is 
inversely proportional to U' and J. Typical values of 
J H eis are 0.14 ±0.01 for (H', J)=(4, 3), 0.042 ±0.002 for 
(V, J)=(16, 13), and 0.019 ± 0.002 for (£/', J)=(32, 28), 
showing that the expected spin-gaps are very small in 
the units that are natural for Hamiltonian Eq.(l). Re- 
ducing J at fixed U', the S^l region is accessed where 
the on-site S=l and S=0 states start mixing, produc- 
ing deviations in the spin correlations from those of the 
Heisenberg model. 

The introduction of mobile holes in the S=l region, 
through the reduction of the electronic density, produces 
the most interesting results of our investigations. In the 
regime of moderate couplings, both U' and J less than 
10, the presence of spin incommensurate tendencies are 
clearly observed (Fig. 2a), in a wide range of densities, 
through the appearance of a prominent peak at a mo- 
mentum different from ir. No fine tuning of parameters 
is needed, the effect exists in a wide range of couplings.lia 
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Note that although experimentsta suggest that the ma- 
terial Y2_ x Ca x BaNi05 remains insulating upcm doping, 
the observed drastic change in the resistivityu justifies 
the use of a metallic description of the compound, as,-pb- 
served in our studies where charge is not localized. □ It 
may occur that holes in the real material are mobile over 
several lattice spacings, but defects destroy the complete 
metallicity of the samples (see also Ito et alSjj). 
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FIG. 3. Spin correlations for one (n=1.9375), two 
(n=1.875), and four (n=1.75) holes, doped on a 16 sites sys- 
tem at n—2. Grey circles indicate the sites where the holes 
were projected at orbital 1. The thickness of the lines is pro- 
portional to the spin bond strengths, according to the scale 
shown. Full (dashed) lines indicate AF (FM) bonds. The cou- 
plings are (U', «7)=(4, 3). In all cases shown, across-the-hole 
bonds are AF and S(k) has spin IC, as in Fig. 2a. 

Insight into the origin of the spin IC can be gained by 
investigating spin correlations in the vicinity of the holes. 
This can be obtained by projecting the hole on a giving 
site out of the ground-state, and then measuring spin 
correlations around and across that hole. This projec- 
tion procedure has been successfully jEtcfiviously applied 
in studies of models for the cuprates.EJEll Fig. 3 shows a 
pictorial representation of a 16 sites system doped with 
one, two, and four holes .□ The holes were projected onto 
their most probable locations on the chain (remember 
that OBC are used). Since hopping is restricted within 
the 3d(3z 2 — r 2 ) orbitals, the projection is always car- 
ried out onto orbital 1. After the projection, the spin 
bond strengths (Sj • Sj) were measured for i and j being 
nearest-neighbor sites, or sites adjacent to the same hole. 
The bond strength is proportional to the thickness of the 
line connecting the sites in Fig. 3. In that figure, the AF 
(Si -Sj) < (FM (Si -Sj) > 0) bonds are represented by 
a full (dashed) line. 

The results in Fig. 3 were obtained for (U 1 , J)=(4, 3), 
but they are representative of tendencies found in most 
cases studied here, for small and moderate couplings. 
Across the projected holes, the clear presence of AF 
bonds is observed. It is this 7r-shift in the staggered spin 
pattern that causes the incommensurability. Intuitively 
these 7r-shifts improve the mobility of holes. Without 
them, the hole movement would lead to the creation of 
"strings" of ferromagnetic bonds, increasing the energy. 
This across-the-hole structure is reminiscent .of, those ob- 
served in models for cuprates in ID and 2DJ3'B with the 
conceptual difference that the present results exist in the 



spin-liquid regime (nonzero spin gap). It appears that 
the local tendency to produce spin IC, in the form of 
across-the-hole AF correlations, is present regardless of 
the spin-liquid vs critical or long-range ordered charac- 
ter of the spin background. Note that in fact for a S=l 
chain it is possible to use the concept of "AF order" in 
the spin background as long as the distances considered 
are smaller than the AF correlation length, which for the 
S=l Heisenberg chain is about six lattice spacings. This 
is sufficiently large to accommodate the small across-the- 
hole structure observed in our studies. This also suggests 
that even if holes are localized in Y2_ x Ca x BaNiQ&_with 
a distorted region of only a few lattice spacings J3'Ej it is 
still possible for holes to generate spin IC. 

The results shown in Fig. 2a are representative of data 
gathered at several couplings, and they not only contain 
information on spin IC, but on other nontrivial spin pat- 
terns as well. In particular, it was observed that for a 
fixed U', the spin IC survives longer the decrease in elec- 
tronic density by reducing J, but eventually a ferromag- 
netic (FM) regime always appears if the hole doping is 
large enough, as shown in Fig. 2a where a prominent peak 
at zero momentum is observed at density 1.5. Then, the 
phase diagram of models for doped S=l chains must con- 
tain a FM phase at low enough electronic density. The 
existence of a competing FM phase is in agreement with 
previous calculations using effective t-J-like models valid 
at low energies, c3 also with expectations coming from 
manganit&irpestigations where the relevant densities are 
close to lJlj'EJ and it could be tested experimentally by 
increasing the hole doping in the S=l chain compounds. 

In addition, doping at stronger couplings than those 
shown in Fig. 2a produces phase separation (PS) between 
hole-rich and hole-poor regions, if J is close to the lower 
boundary of the S=l region of Fig. la. This is the case 
for the results shown in Fig. 2b. Such PS regime is char- 
acterized by a negative inverse compressibility, the pres- 
ence of two simultaneous peaks in the spin structure fac- 
tor at a fixed total electronic density, and a tendency to 
charge clustering in the local density. Tendencies toward 
FM and electronic PS found here are similar to those 
encountered in Kondo-like models for the manganites,Eil 
and effective models for nickelates.E^I Manganite models 
also have two orbitals, but in addition they contain local- 
ized tig spins. In the S=l chains considered here, the or- 
bital (x 2 — y 2 ) is fully occupied, and plays a role analog to 
that of those localized spins of Mn-oxides. Then, the ten- 
dencies to ferromagnetism observed here may be caused 
by a mechanism similar to Double Exchange, which could 
be operative in highly-doped Ni-oxides. 

The analysis of the previous paragraphs leads to a sec- 
ond possible rationalization of the spin IC. As explained 
before, hole doping eventually drives the system into a 
FM phase. But the transition between the IC phase dis- 
cussed previously and the FM phase is not always abrupt. 
Fig. 4 contains results analogous to those of Fig. 3 but 
now for ([/', J)=(8, 6). In all three doped systems dis-. 
cussed, S(k) was still found to be incommensurate.EHl 
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In Fig. 4, the spin correlation pattern when one or two 
holes are projected to their most likely position is sim- 
ilar to Fig. 3, but for the case of four holes the spin 
configuration is clearly different. Instead of AF bonds 
across the holes, a small three-site island of ferromag- 
netism is formed around each hole, resembling a mag- 
netic polaron. The presence of FM polarons was already 
observed in recent studies of doped two-orbital chains 
with strong Hund's rule couplings l3 These polarons are 
coupled in such a way that they generate an overall spin 
IC (in particular note the way in which the two central 
polarons are coupled antiferromagnetically) . Then, there 
is a second mechanism to generate spin IC in this con- 
text. However, since experimentally no evidence of FM 
has been reportcd,B the first explanation with AF corre- 
lations across holes appears more realistic. 
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FIG. 4. Same as Fig. 3 but for (U', J)=(8, 6). Although 
in all cases shown S(k) is still incommensurate, the system 
doped with four holes resembles a polaronic lattice, with the 
holes located at the centers of small FM islands. 

Note that the presence of phase-separation tendencies 
reveal a source ofj-attraction among holes. Actually, in 
previous literatures it was already shown that holes at- 
tract in a large region of parameter space, using t-J-like 
models for Ni-chains. This is reasonable, since each hole 
creates a spin distortion, and sharing distortions reduces 
the energy. The pairing of holes in a spin-gapped sysfc. 
tern has before been documented in models for ladders, — 
and it appears in some models for doped S=l chains!^ 
Studies of two-orbital models have also revealed robust 
singlet-pair correlations.^ Then, it is tempting to spec- 
ulate that the material Y2- x Ca x BaNi05 may become 
superconducting upon the application of pressure, as it 
happens with the two-leg ladder compound [14-24-41]. 
Experimental studies of doped S=l chains under high 
pressure may lead to interesting results. 

Summarizing, studies of hole-doped S=l chains show 
that spin incommenpirate correlations can be generated 
in a realistic model,E-3 as observed in experiments. The 
mechanism is similar to the phenomenon recently dis- 
cussed in S=l/2 ladders and planes, namely the genera- 
tion of AFpCorrelations across holes to optimize the hole 
movement!] It is expectedlhat this mechanism may con- 
tribute, at least in part,Ot3 to the understanding of the 
recently unveiled spin IC tendencies in Y2- x Ca x BaNiOs. 
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